Fermi level alignment in molecular nanojunctions and its relation to charge transfer 
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The alignment of the Fermi level of a metal electrode within the gap of the highest occupied 
(HOMO) and lowest unoccupied orbital (LUMO) of a molecule is a key quantity in molecular 
electronics, which can vary the electron transparency of a single molecule junction by orders of 
magnitude. We present a quantitative analysis of the relation between this level alignment (which 
can be estimated from charging free molecules) and charge transfer for bipyridine and biphenyl 
dithiolate (BPDT) molecules attached to gold leads based on density functional theory calculations. 
For both systems the charge distribution is defined by a balance between Pauli repulsion with 
subsequent electrostatic screening and the filling of the LUMO, where bipyridine loses electrons to 
the leads and BPDT gains electrons. As a direct consequence the Fermi level of the metal is found 
close to the LUMO for bipyridine and close to the HOMO for BPDT. 

PACS numbers: 73.63.Rt, 73.20.Hb, 73.40.Gk 



Interest for electron transport in nano-scale con- 
tacts has recently intensified, because (i) the advent 
of the technologically motivated field of molecular elec- 
tronics^'^'^'^ (ii) recent progress in the experimental 
techniques for manipulating and contacting individual 
molecules^&L^, and (iii) the availability of first princi- 
ples methods to describe the electrical properties of sin- 
gle molecule junctions. ^'^^'^^'^^'^^ These latter methods 
are usually based on density functional theory (DFT) 
in combination with a non-equilibrium Green's function 
formalism.^^. 

It has become clear through many studies that the 
conductances of molecular junctions can be entirely con- 
trolled by the positions of individual molecular levels. 
In for example the case of a molecular contact consist- 
ing of a bipyridine molecule attached to Au leads it has 
been demonstrated^^' that the transmission is com- 
pletely dominated by the position of the lowest unoc- 
cupied molecular orbital (LUMO) and that this position 
may vary significantly with for example the surface struc- 
ture of the leads. In order to describe the transport in 
molecular junctions it thus becomes a key issue to deter- 
mine factors controlling the line-up of molecular levels 
relative to the Fermi level of the metal. Is the alignment 
governed mostly by the decay of the potential at the sur- 
face or it is controlled by level shifts due to the specific 
interaction between a particular molecule and a particu- 
lar surface? In order to approach this question we study 
in the following two different molecular junctions which 
turn out to behave quite differently with respect to level 
alignment. 

Fig. n compares bipyridine and biphenyl dithiolate 
(BPDT) attached to gold electrodes with the same sur- 
face structure and the same bonding configuration (top 
panels). In both cases no molecular levels can be found 
close to the Fermi energy if the vacuum potentials from 
separate calculations for the molecules and the metal sur- 
face are set equal (middle lower panel). However, when 
the molecules are coupled to the electrodes, there are 



peaks in the transmission functions (left and right lower 
panels), where Ep is crossed by their tails at different 
sides of a gap for the two different molecules. For bypri- 
dine, the LUMO has moved downwards with respect to 
its original (non-interacting) position, for BPDT it has 
moved upwards. The exact position of these peaks de- 
termines the conductance which is defined by the value 
of the transmission function at Ep and can vary between 
0.03 and 0.44 Go for BPDT depending on the bonding 
configuration. What we shall show in the following is 
that the level shifts can be directly determined from an 
appropriately defined charge transfer between the metal 
surface and the molecule. The down-shift of the bipyri- 
dine LUMO is for example associated with an electronic 
charge transfer from the molecule to the surface. This 
might at first seem counterintuitive since the down-shift 
of the LUMO apparently leads to a slight occupation of 
this state (see Fig. but still the net electron transfer 
is away from the molecule due to a combination of Pauli 
repulsion and screening effects as we shall demonstrate. 
For benzene dithiolate (BDT) and BPDT the alignment 
of molecular levels coupled to a Au (111) surface has 
previously been suggested to be linked to charge being 
transferred from the surfaces to the molecules. 

All electronic structure calculations in this study are 
performed using a plane wave implementation of DFT— 
with an energy cutoff of 340 eV, where we used ultra- 
soft pseudopotentials^^, and a PW91 parametrization for 
the exchange and correlation functional. The transmis- 
sion functions of the molecular junctions in Fig. Q were 
calculated using a general non-equilibrium Green's func- 
tion formalism for phase-coherent electron transpor1i2^, 
where both, the Green's function of the scattering re- 
gion and the self-energies describing the coupling to the 
semi-infinite electrodes, were evaluated in terms of a ba- 
sis consisting of maximally localized Wannier functions 
obtained from a transformation of the Kohn-Sham eigen- 
states^^ within an energy range up to 2 eV above Ep. In 
our calculations the supercells for the scattering region 



2 




FIG. 1: (Color online) A bipyridine (left) and BPDT (right) molecule suspended between Au electrodes. The middle lower 
panel shows the alignment of the molecular levels with the metals Fermi energy by equalizing vacuum potentials of the isolated 
molecules and surfaces where Ep, and the HOMO and LUMO are marked with gray (red) lines. The left and right lower panels 
show transmission functions for the coupled systems for bipyridine and BPDT, respectively. 



are defined by 3 x 3 atoms in the directions perpendic- 
ular to the transport direction and contain three to four 
surface layers on each side of the molecule. We found 
that a 4 X 4 k-point grid is needed for the sampling in 
the transverse Brillouin plane in order to obtain well con- 
verged results for the conduct ance^i. 

We investigated the variation of the energetic position 
of the lowest-lying molecular orbital (denoted MOl in the 
following text) with respect to the metal's Fermi level, in 
dependence on the distance between the surface and the 
molecule d. This is depicted in Fig. [SJ where do marks 
the equilibrium bond length between the nitrogen atoms 
of the bipyridine molecule (Fig. ^) or the sulfur atoms 
of BPDT (Fig. Eb) and the Au atoms they are attached 
to. Since MOl is ~ 10 eV below the lowest-lying Au 
valence states, its energetic position must be exclusively 
guided by rigid potential shifts without any direct hy- 
bridisation effects. Now we want to address the question 
whether these rigid potential shifts have a quantifiable re- 
lation to charge being transferred between the molecule 
and the surface. In our study we make use of the con- 
cept of fractional charges, which can be introduced either 
through a straight forward extension of the mathemat- 
ical framework of DFT or via a thermodynamic ensem- 
ble interpretation of the mixing of pure states. This 
makes it possible to determine the ground state electron 
density and electronic eigenvalue spectrum for a (albeit 
only finite) system with fractions of electrons removed or 
added when compared to the total charge of all the nu- 
cleis. The shift of MOl in the charged isolated molecules 
is shown as insets in Fig. [21 for bipyridine and BPDT. 
As can be seen the molecular levels move up when the 
molecule is charged as a consequence of electrostatic re- 



pulsion. Our main argument is that these level shifts of 
the free molecule can explain the level shifts in the cou- 
pled molecule surface system. To show this we use the 
charge density differences between the coupled and iso- 
lated subsystems and the linear energy /charge relation 
for the free molecule. This allows us to calculate the po- 
sition of MOl as a function of distance, which we then 
compare directly with the actual energies from the cou- 
pled system in Fig. [3 For a physical interpretation of 
such fractional charges two notes of caution have to be 
made: i) Fractional charges are only physical meaningful 
for small distances d, where an equilibrium between the 
molecule and the surface can be assumed. Our plots for 
long d are only used for simulating a gradual switching on 
of the interaction, ii) In our calculations we do not cor- 
rect for the lack of the derivative discontinuity (DD) of 
the exchange correlation potential within standard DFT 
methods, which can lead to an underestimation of the 
HOMO-LUMO gap and an overestimation in the conduc- 
tance for weakly coupled single molecule junctions even 
at low biases^^. Including DD would make our scheme 
impractical, since MOl would jump discontinously when 
going from negative to positive charges. However, our 
study does not focus on absolute values for conductance 
or gap size but on the comparison of equilibrium charge 
transfer in two different junctions. 

For bipyridine at large distances MOl rests at an 
energetic position eyac-aiigm which corresponds to the 
one it would hold if the vacuum levels of the isolated 
molecule and surface were aligned (Fig. [^J. This indi- 
cates that there is no interaction between the two sub- 
systems and that both are entirely charge neutral for d 
larger than ~ 6 A . At the bonding distance do an ef- 
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FIG. 2: (Color online) MOl energies (relative to Ef of the 
Au surface) depending on the distance d between the molecule 
and the Au surfaces for a) bipyridine and b) BPDT. The black 
solid curves are taken from the coupled systems, for the red 
curves estimates from the charged free molecules have been 
used, and the dashed line shows the position from vacuum 
alignment without charging (evac-aiign)- MOl in its depen- 
dence on the charge for the free molecules is shown as insets. 



fective charge of -0.25 electrons on the bipyridine can be 
derived from the shift of MOl by comparison with the 
charged free molecule. Partial charges with respect to the 
isolated subsystems can also be computed directly from 
electron density differences, which results in a charge of 
-0.23 electrons on the bipyridine molecule at the same 
distance. As can be seen from Fig.Q] the HOMO-LUMO 
gap is much smaller for BPDT (~ 0.3 eV) than for bipyri- 
dine (~ 3.2 eV). This is because BPDT lacks two elec- 
trons, which are subtracted from its tt system, when two 
hydrogen atoms are removed from the stable aromatic 
molecule biphenyl dithiol in order to form the highly re- 
active biradical BPDT, which is then attached to the Au 
surfaces. Since the molecular levels corresponding to the 
dangling bonds on the sulfur atoms are fully occupied in 
BPDT, the HOMO of biphenyl dithiol gets emptied and 
becomes the LUMO of BPDT. A further difference be- 
tween bipyridine and BPDT (see Fig. is that the long 
distance position of MOl for BPDT is not just eyac-aiign 
but is instead ~ 0.5 eV higher in energy. Vacuum level 
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FIG. 3: (Color online) Charge density difference An{xz) 
(black solid line, summed up parallel to the surface plane) 
and its integral (black dashed line) for bipyridine coupled to 
Au atoms ((i=2.63 A ). An{xz) constructed only from contri- 
butions of the bipyridine HOMO (see inset) and the Au s and 
d (z^) states is also given for comparison (green line). 



alignment for BPDT on Au (111) leads to a situation 
where the LUMO lies well below Ep (see Fig. QJ, and 
independent of d the two subsystems cannot exist within 
the same cell without charge being transferred from the 
surface to the molecule. We find 0.11 electrons gain on 
the BPDT molecule from the MOl level shifts and 0.09 
from charge density differences An{xz) at large d. At 
the bonding distance do, the BPDT molecule gains 0.12 
electrons calculated from MOl level shifts and An{xz) 
cannot be interpreted unambigously due to the strong 
hybridisation of Au and S states. We note that BPDT 
differs from bipyridine in that charge transfer occurs in 
the opposite direction, where the amount of charge being 
moved is by almost a factor of two smaller for BPDT and 
its dependence on d is negligible. 

An apparent problem with the picture we propose here 
is the following: In the case of bipyridine charge is mov- 
ing away from the molecule and the levels are therefore 
shifting down in energy. However, as the LUMO hits 
Ei? it must begin to fill and this is naturally associated 
with charge being transferred to the molecule. Both ef- 
fects are in fact taking place. For the further analysis 
of the bipyridine junction, we use a model system where 
the leads are replaced by single Au atoms. Looking at 
An{xz) in Fig. |31we find maxima at the Au atoms and 
minima close to but not at the nitrogen positions. When 
we form charge density differences just using the bipyri- 
dine HOMO (which is twice degenerate and both orbitals 
are fully occupied), the Au s and the also fully occupied 
d (z^) states, the result still has the same nodal structure 
as An{xz). The two functions differ in the minimum at 
the nitrogens being deeper for the difference constructed 
only from six orbitals. In the latter case there is also 
no charge depletion in the center of the molecule. Both 
aspects can be explained in terms of screening. Since a 
very localised minimum in the charge density is energet- 
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ically unfavourable, the lower lying MOs are polarised 
so that this minimum is at least partially smoothed out. 
Such a polarisation of MOs has the effect that the total 
charge density at the center of the molecule is reduced, 
which explains the net charge transfer from the molecule 
to the Au. For the maximum at the Au atoms no such 
screening occurs since for the bare atoms there are no 
electrons available for achieving that. For the realistic 
surface calculations, however, the situation is different 
and the same peak in An{xz) (not shown here) becomes 
smaller and broader. 

But what is the effect of the partial filling of the 
LUMO, which occurs in the bipyridine junction only for 
small d? It can be seen in the change of curvature in the 
MOl energies vs. d in Fig.[2jL for distances smaller than 
di. Because now the process described in the last para- 
graph is partially neutralised by charge flowing back to 
the molecule, the lowering in energy of MOl is slightly 
reduced, thereby flattening the curves. For Fig. [31 we 
have chosen a distance larger than di in order to be able 
to discuss the two driving forces for charge movement 
separately. For BPDT the same balance between Pauli 
repulsion of electrons (leading to electron depletion on 
the molecule) and the filling of the LUMO exists but the 
latter process dominates, and the net effect is therefore a 
movement of electrons towards the molecule. The strong 
system dependence of the Fermi level alignment in molec- 
ular junctions can also be seen by a comparison of our 
results with a recent article on alkane thiolates on Au 
(111)^^, where the dominant effect was found to be due 
to molecular dipoles which result from the asymmetry of 
the system. 

In summary, we presented a detailed analysis of the 



energetic alignment of molecular orbitals with respect 
to the Fermi level of gold electrodes in single molecule 
nano junctions with bipyridine and biphenyl dithiolate 
molecules. The outcome of this alignment has a cru- 
cial effect on the zero bias conductance of the junction, 
the major source determining it is equilibrium charge 
transfer between the molecule and the electrodes. We 
established that a comparison between the energies of 
the lowest lying molecular orbitals within the junction 
and for the isolated but partially charged molecule gives 
reasonable estimates for the net charge being transferred. 
For the bipyridine, the molecule is drained of electrons, 
the biphenyl dithiolate on the other hand gains electrons 
from the Au electrodes. From a simplified model of the 
junctions where the electrodes have been replaced by sin- 
gle Au atoms, we could derive that the charge transfer 
process for the two molecules we studied is determined 
by the balance of two effects. Pauli repulsion between oc- 
cupied molecular and Au orbitals pushes electrons away 
from the molecule and the interface region to more re- 
mote parts of the electrodes surface or bulk, which fol- 
lowed by screening depletes the molecules of electronic 
charge. The filling of the LUMO on the other hand re- 
sults in an electron surplus on the molecule. For bipyri- 
dine the first effect dominates, for BPDT the situation is 
reversed. 
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